Distribution of Transmitted Charge through an Ultrasmall Double- Tunnel Junction 
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The transmission of charge through an ultrasmall double junction is considered with Coulomb 
effects but at zero temperature. We construct an equation which describes the time development of 
the transmission probability of charges and solve this equation in terms of a recursion relation. The 
results are compared with the double junction without charging effects. 

73.23.Hk, 73.40Gk, 73.40Rw 



I. INTRODUCTION 



II. TIME DEVELOPMENT 
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During more than one decadeH'tJ single-electron tunnel- 
ing (SET) has experienced considerable progress in both 
experiment and theory. Basically, SET is a stochastic 
tunneling process of charges occurring at ultrasmall low- 
capacitance high-resistance tunnel junctionsfl'Q Those 
tunnel junctions ensure that, on the one hand, the charg- 
ing energy is large (compared to the thermal energy) 
whereas, on the other hand, charge fluctuations are sup- 
pressed. Thus, the number of tunneled charges forms 
a classical variable whereas the tunneling probabilities 
(rates) are calculated quantum mechanically, for instance 
in terms of Fermi's Golden Rule.H 

This calculation restricts the application of the semi- 
classical approach to slow time dependencies on the scale 
of the iiC-time of the tunnel junction. It is in the range 
of 10~ 10 s for nowadays devices. There are approaches be- 
yond the semi-classical theory that allow the exploration 
of short time scales (for instance Ref. ^3j) , however for the 
purpose of this work simpler Golden Rule approach is 
sufficient. 

This paper is devoted to the study of the ultrasmall 
double tunnel junction, which is one of the cornerstones 
of single electronics in experiment and theory. It con- 
sists of two junction connected in series, thus forming a 
small island inbetween. This system is capable of display- 
ing pronounced effects of single-electron tunneling like 
Coulomb blockade, Coulomb staircase, and shot-noise 
suppression for suitable parameters of the applied bias 
voltage V. 

Throughout this paper zero temperature is assumed. 
On the one hand, this assumption simplifies our consider- 
ation considerably and enhances its clarity. On the other 
hand, many SET experiments are done at sub-Kelvin 
temperatures, where a zero-temperature approximation 
is reasonable. 

In the next Section we develop the basics of our ap- 
proach and display the connection to the well-known "or- 
thodox" theory and to the double junction without charg- 
ing effects. It follows a discussion of these results in the 
third Section. 



For the double junction 
suits in a master equation 



semi-classical approach re- 
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dt 



r\ (m — l)<r(m — 1, t) 

+r 2 (m + l)a(m+ l,t) 
— [ri(m) + r 2 (m)]cr(m,i), 
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which describes the time dependence of the probability 
a(m) of m extra charges on the island. r\$ describe the 
rates of a charge tunneling event at the first and second 
junction. Tunneling against the bias voltage does not 
occur for zero temperature. The number of the charge 
states m is finite, their energy lies between the Fermi 
levels of the banks. We denote the lowest and highest of 
these states by min and max, respectively. 

With regard to the stochastic process of single-electron 
tunneling, the first moment describes the average num- 
ber of tunneled charges, (n(t)}, which is connected with 
the stationary electrical current through the double junc- 
tion, (I) = e(n(t))/t, whjere e_is the unit charge. Using 
the stationary solutiorM3 o 



a(m) 



of ® 
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with the normalization constant Z, 
rent reads 



the stationary cur- 



(7) = e r\(m)cr(m) = e r2(m)a{m). (3) 



m— mm 
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We note that continuity levels the current across the first 
and the second junction. 

The second moment of the stochastic process under 
consideration is connected with the current noise in a 
double junction. |— It has been investigated theoretically 
in recent yearsETEJJ and its low-frequency limit was mea- 
sured in experiment. Ej 

Further moments were - to our knowledge - not sub- 
ject of theoretical efforts, however, an approach to the 
stochastic process of tunneling through a double junc- 
tion (without charging effects and at zero temperature) 
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was performed recently.lll It is based on the calculation 
of the probability P n (t) of n transmitted charges in terms 
of the equations 



dP (t) 

dt 
dP n (t) 

dt 



= -72^0 (t) 



(4) 



71.2 are the transmission probabilities of the two junc- 
tions. According to the Pauli exclusion principle there 
are only two different island states that have to be con- 
sidered, unoccupied (0) and occupied (1) by one charge. 
The initial occupation probability of them is given bylla 



P 



(0,1) 



(t = 0) = 



72,1 
7i + 72 



and the probability of n transmitted charges at arbitrary 
time t i£3 



Pn(t) = 



72 



71 +72 



7i 



7i +72 



This problem is solved analytically in Ref. [T^, i.e. the 
transmission probabilities P n {t) can be calculated. The 
description of the stochastic process, however, is ipych 
clearer in terms of its characteristic function x( A, t)L2l 



P n {t) e 



i n A 



(5) 



n=0 



from which in turn all moments /ifc(t) = (n(t) k ) of the 
stochastic process are obtained 

Equivalently, the cumulants follow from the logarithmic 
expansion of the characteristic function x(A,i).E3 

In our approach we include charging effects, however, 
still assume zero temperature. Therefore, a finite number 
of island charge states m is considered and the descrip- 
tion is done in terms of the probability p(m, n, t) of n 
transmitted charges through both junctions and simul- 
taneously to excess charges on the island between them. 
Thus, the following sum rules hold 
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1= ^2 X^( m,n '*)' 

rn— min n— 

oo 

a(m,t) = y^p(m,7i,t), 

n=0 
max 

Pn{t) = Y p(m,n,t). 



(0) 
(7) 
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For the time development of p(m, n, t) the following 
first-order differential equation is found 



dp(m, n, t) 
dt 



= ri(m — l)p(m — 1, n, t) (9) 

+r2{m + l)p(m + l,n — l,t) 
— [ri (to) + r 2 (to)] p(m, n,t). 



The discussion of this equation and its solution is the 
main point of this paper. 

The initial condition for the solution of (^|) follows from 
(0) using P (t = 0) = 1, 



p(m, 0, t = 0) = £t(to) 
p(to, n, t = 0) = 



(10) 



n > 0. 



Summing (|^) over n and m, respectively, results by 
the use of (0) and (||) in the former Eqs. [l] and |[ Even 
(^) includes Coulomb effects in our case indicated by to 
entering into (||). 

Let us consider ([)]) in the long-time limit (t — > oo) 
shortly. Single-electron tunneling is usually considered 
in this limit and several simplifications apply. The prob- 
ability p(to, n, t) factorizes p(m, n, t) — P n (t)o(m) 1 where 
aim) is the stationary solution (||) of the master equation 
(fy). Furthermore, the transmission probabilities become 
periodic, P n±1 (t) = P n (t =F e/(/)). Under this conditions 
the stationary-current formula ([5]) results from (^|), even 
for non-zero temperatures. 

In general, the assumptions of the last paragraph do 
not hold and the solution takes a more complicated form, 
however, it can be expressed semi-analytically in terms 
of a recursion formula. The general solution of (|^) can 
be written as 



p{m,n,t)= J2 A l k m-^M-lkt) (11) 

A:— min 

n—1 max 

+ J2 ^nm7fexp(-7 fe t), 



/— A;— min 



where 7^ = ri(k) + r2(k) is used for abbreviation. Thus, 
the solution turns out to be a superposition of different 
poisson-like probabilities. The stochastic process as a 
whole, however, is not poisson-like. For n = the second 
term in (|ll]) vanishes. 

The coefficients A l £ m have to fulfill a recursion relation 
which can be derived from (||) by introducing (11). For 
a short-hand notation of this recursion relation we make 
use of the following decomposition 



1 = gj+lrn 

(S + 7m) (« + 7fc)' +1 S + 7m 



Ik 



which defines the numbers a\, n in a unique way as solu- 
tion of a linear equation system. In terms of this notation 
we find 
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Kt = ri(m- mH^alt k < m (12) 

KZ = n{m- l^- 1 - + r 2 (m + 

m— 1 n — 1 max 

^ = n(m-l) £ ^-i<li m + n(m-l)^ £ 4L-i<i m 

fc=mill 1 — fe = min 

m+l n— 2 max 

+r 2 (m + l) £ ^:}" l+ l< m lfc +r 2 (m + l)^ A n-lm+M k +lm 

fc = mm ; = *=mla 

Aun =n{m- l)A 1 -^ + r 2 (m + lji^+i 1 < i < n - 1 

n-l 

<™ = ri(m - l^.^If [m > fc] + n(m - 1) £ < fc ro _ lC *^ 

+r 2 (m + l)A» = ^ +1 <- lfe If [™I > ] +r 2 (m + 1) £ A* lm+1 <. 



For simplicity we made use of the Iff. . .] notation: the 
corresponding terms contribute only if the given condi- 
tion is fulfilled. Eq. |lj is, in connection with (|ll|), our 
main result. 

In addition to (^2|), the initialization of A l £ m follows 
from (||) for n = using the given initial condition (|io|), 



A 



Ok 



r x (m-l)>C-i 



^ = < 7(m) + r 1 (m-l) X 1 



fc < m (13) 



, * Ik ~ 7m 
fc— min 

Eq. [l^ does not hold in the very symmetric case. Then, 
the solution might derived from the given set of equation 
in terms of a limit procedure 7^ — > 7 m . 

Eqs. [l^ and |l^ allow the construction of all A l £ m - 
These coefficients tell the important contributions within 
( |ll| ) from the unimportant ones. If a sole term dominates 
the stochastic process will approach a Poiaaon process, 
similar to the case without Coulomb effectsJl3 In general, 
however, more terms contribute and the total process is 
a superposition of Poisson processes. The result (12), 
however, does not seem to provide clear guidance to a 
simplification by the selection of relevant terms. Even if 
the given recursion relation (|l^) looks clumsy its numeric 
implementation is neat. 



III. DISCUSSION 

Owing to the complicated expressions ( |l2| ) we discuss 
our results by means of numerical computations. The 
most important question concerns the time domain where 
( |TT| ) is valid. Due to the derivation of the rates ri, 2 (m) 
by Fermi's Golden Rule, the time scale is limited to 
t ^> i?i, 2 Ci i2 . On the other hand, for very long times 
the system will approach the mentioned stationary state 
of the island and the investigation looses its thrill. 
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FIG. 1. Transmission probability of charge through a dou- 
ble junction with (solid lines) and without (dashed lines) 
Coulomb effects for the first ten transmitted charges. The 
inset shows displays the sum of these ten probabilities. 

In Fig. [l] we show the results for P n {t) for n = ... 9. 
The parameters of the double junction under investiga- 
tion are R x = R 2 = lOOkft, C x = lOaF, C 2 = lOlaF. The 
solid lines show the data corresponding to (0), whereas 
the dashed lines reflect the results without consideration 
of the charging effects according to Ref. [I4 

Since not shown in the figure let us mention the ini- 
tial condition P (t = 0) = 1, which is fulfilled by both 
sets of transmission probabilities. The suppression of the 
charge transmission by Coulomb interaction, however, is 
clearly observed in Fig. |[ The probability P n (t) with 
charging effects extends over a wider time range than 
the same probability without charging for a given num- 
ber n of transmitted charges. Thus, for a fixed time t 
less charges are transmitted in former case and the con- 
duction is suppressed. 

In the inset the probability 
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71=0 

is shown. For short times it approximates the total prob- 
ability 1, whereas later neglected states get more impor- 
tant. We note that the first ten P n (t) are sufficient to de- 
scribe times up to « 20 Ry,Cy, in the case with Coulomb 
interaction, but only up to « 2 i?sCs in the case without 
charging. For short times the applied description of the 
system fails, however. 
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FIG. 2. Fano factor r(t) for a double junction with and 
without Coulomb interaction (solid and dashed line, respec- 
tively). The short-dashed line corresponds to the long-time 
limit of the case without charging effects. 



In order to apply the solution 
factor 



Llh , we plot the Fano 



r(t) 



for the case with and without Coulomb interaction in 
Fig. ||. The Fano factor describes the correlation of sin- 
gle tunneling events. Small r(t) indicate a large degree 
of correlation. Therefore, a high Fano factor is expected 
at t = 0, whccnas lateron it should approach a lower sta- 
tionary limit .Era 

Again, our computation makes use of the first ten P n (t) 
only. This corresponds to Fig. Ill Owing to the limited 
number n the plotted approximation fails at longer times, 
i.e. when the total probability (as shown in the inset of 
Fig. [I]) differs from unity. Besides this flaw the expected 
behavior of the Fano factor is recovered. For the case 
without Coulomb interaction we plot the long-time limit 
of Ref. [l| as well. 

Coulomb interaction favors correlated tunneling. This 
is seen by the lower value of r(i) in comparison to the sta- 
tionary case without Coulomb interaction. The degree of 
correlation depends on the symmetry of the double junc- 
tion [...]. It gets maximal in case of a symmetric double 
junction (i?i = i?2, C\ = C%), where the Fano factor 
reaches 1/2. The double junction under consideration, 



however, is rather asymmetric. Thus, the decrease of the 
Fano factor is small. 

For the studied case with Coulomb interaction notice- 
able changes in r(t) are observed even in the time range 
of t k 10-RsCs. This is well above the critical short- 
time range, where the approach in terms of the Golden 
Rule is insufficient. This fact supports justification of the 
presented approach. 
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FIG. 3. Absolute value of the characteristic function 
|x(A,t)| for three different values of A (0.5, 1.0, 2.0, from top 
to bottom). The double junction parameters correspond to 
the other figures. Again, solid and dashed lines indicate the 
computation with and without Coulomb energy, respectively. 

Finally, in Fig. || we display the time-dependence of the 
characteristic function %(A, t) of Eq. |[ As mentioned 
above, it "characterizes" the stochastic process of elec- 
tron tunneling through the double junction. Again, we 
use the probabilities of the first ten transmitted charges 
for computation. Therefore, the relevant time range is 
here similar to that of Fig. |l| and Fig. |^. This plot shows 
that the presented recursion relation scheme can be used 
to compute further properties of the underlying stochas- 
tic process. 



IV. CONCLUSION 



In this paper we have studied the stochastic process of 
single-electron tunneling in an ultrasmall normal-metal 
double-tunnel junction at zero temperature. We have 
introduced an equation, (^), which enables the investiga- 
tion of transmission probabilities of charge correspond- 
ing to the "orthodox theory" of single-electron tunneling. 
The time-dependent solution of this equation is derived in 
terms of a recursion relation, (|l^) , with the appropriate 
initialization The results are discussed in compar- 

ison to the known investigation of the double junction 
without charging effects.Ej We observe a suppression of 
tunneling due to Coulomb repulsion and an enhanced 
correlation. The recursion relation provides a suitable 
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tool for computation of stochastic properties of the dou- 
ble junction. The absolute value of the characteristic 
function x(A, t) serves as an example. 
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